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Abstract 

In  a  high  average  power  free  electron  laser  (FEL)  the  wall-plug  efficiency  is  of  critical 
importance  in  determining  the  size,  complexity  and  cost  of  the  overall  system.  The  wall- 
plug  efficiency  for  the  FEL  oscillator  and  amplifier  (uniform  and  tapered  wiggler)  are 
strongly  dependent  on  the  energy  recovery  process.  A  theoretical  model  for  electron 
beam  dynamics  in  the  energy  recovery  linac  is  derived  and  applied  to  the  acceleration  and 
de-acceleration  of  nano-Coulomb  electron  bunches  for  a  tapered  FEL  amplifier.  For  the 
tapered  amplifier,  the  spent  electron  beam  exiting  the  wiggler  consists  of  trapped  and 
untrapped  electrons.  De-accelerating  these  two  populations  using  different  phases  of  the 
radio-frequency  wave  in  the  recovery  process  enhances  wall-plug  efficiency.  For  the 
parameters  considered  here,  the  wall-plug  efficiency  for  the  tapered  amplifier  can  be  ~ 
40%  using  this  method. 
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I.  Introduction 

The  purpose  of  this  paper  is  to  obtain  estimates  of  the  wall-plug  efficiencies  for 
high-power  free  electron  lasers  (FELs)  employing  energy  recovery  hnacs  (ERL).  The 
FEL  wall-plug  efficiency  can  be  substantially  increased  by  use  of  ERLs  [1].  For  high- 
power  FELs,  the  wall-plug  efficiency  is  of  major  importance  in  determining  the  overall 
size  and  cost  of  the  system.  The  power  flow  diagram  for  an  FEL  based  on  an  ERL  is 
shown  schematically  in  Fig.  1 .  In  Fig.  1  the  injected  electron  beam  is  accelerated  by 
radio-frequency  (RF)  fields  in  the  ERL  and  sent  to  the  FEL  interaction  region.  The  spent 
electron  beam  from  the  FEL  is  sent  back  to  the  RF  section  and  decelerated,  thereby 
replenishing  the  RF  field.  The  decelerated  beam  is  then  sent  to  a  beam  dump.  In  what 
follows,  we  define  the  wall-plug  efficiency  and  calculate  it  for  the  FEL  oscillator  and 
amplifier  (uniform  and  tapered  wiggler)  configurations. 

II.  Wall-Plug  Efficiency 

In  the  steady  state,  the  following  power  balance  relationship  holds  for  the  various 
powers  flowing  through  the  closed  dashed  curve  in  Fig.l, 

n (  Pwp  -  Pcrw  )  =  Pfel  +  Pdump  .  where  V  =  0  -  +  £71rp  is  a  combined  efficiency, 

Pwp  is  the  wall-plug  power,  P(Vro  is  the  power  for  the  cryogenics,  magnets,  etc.,  PFEL  is 
the  FEL  radiation  power  and  Pdump  is  the  electron  beam  power  entering  the  dump.  The 
combined  efficiency  rj  represents  the  efficiency  of  generating  the  RF  fields  for  the 
injector  and  ERL  as  well  as  for  generating  the  injected  electron  beam.  Here  T}rf  is  the 
efficiency  of  generating  RF,  e  is  the  fraction  of  power  going  into  RF  generation  and  rj 
is  the  efficiency  of  the  injector  (see  Fig.  1).  Note  that  the  efficiency  terms  (1  -  e )rjmj  and 
stjrf  are  not  independent  of  each  other.  The  acceleration  factor  in  ERL  is  given  by 
1  +  £J]rf  /(I  -  £)Tj ,  which  is  assumed  to  be  given.  For  example,  if  the  initial  electron 
beam  is  accelerated  by  a  factor  of  15 ,  then  the  ratio  £J]rf  /(I  -  £)t)m]  =14 .  Using  the 
above  power  balance  relationship  the  wall-plug  efficiency  is  given  by 

hwP  =  ^FEL^^WP  =  ^^FEL  »  (1*0 
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where  the  gain  factor  is 

G  =  - - - ,  (lb) 

HfEL  ^dump !  V  ^cryo  I 

Ph  is  the  electron  beam  power  entering  the  FEL  wiggler  and  77fel  =  Pm  /  Ph  is  the  FEL 
interaction  efficiency.  In  obtaining  the  wall-plug  efficiency,  synchrotron  radiation  losses 
and  halo  formation,  among  other  things,  have  been  neglected.  These  effects  can  be  easily 
included  when  estimating  the  system  efficiencies.  Typical  values  for  the  combined 
efficiency  are  r)  ~  0.4  —  0.6  and  77  Pcrm  /  Ph  ~  10“2 .  In  this  paper  we  take  the  average 

beam  current  to  be  Ih ,  the  beam  energy  Eh  =  80MeV ,  Emm  —  3MeV ,  and  Pcryo  =  1  MW 
( for  lh  =  1 A ).  In  general,  the  wall-plug  efficiency  can  be  far  higher  than  the  FEL 
interaction  efficiency,  J]FEL ,  when  an  ERL  is  employed,  i.e.,  G  >  1 .  The  wall-plug 
efficiency  however  is  always  less  than  77 . 

The  beam  power  entering  the  beam  dump,  Pdump ,  depends  on  the  FEL  interaction, 

phase  space  manipulation  of  the  electrons  upon  exiting  the  wiggler  and  the  subsequence 
deceleration  in  the  ERL.  The  deceleration  in  the  ERL  depends  on  the  RF  phase,  self 
fields,  longitudinal  and  transverse  emittances  as  well  as  external  focusing  fields.  These 
effects  will  be  address  in  detail  in  Sections  III  and  IV. 

Within  the  FEL  interaction  region  the  electron  beam  loses  energy  and  develops  an 
energy  spread.  In  addition  to  this  induced  energy  spread  the  beam  has  an  initial  intrinsic 
energy  spread  which  is  typically  <  1  %  of  the  beam  energy  and  is  neglected.  There  is  a 
minimum  energy,  Emin ,  that  the  electrons  can  be  decelerated  to  in  the  ERL.  This 
minimum  energy  limit  is  due  to  various  effects,  such  as,  excessive  electron  de-phasing  in 
the  RF  field  and  transverse  emittance  blowup, 
a)  Oscillator 

Consider  first  an  untapered  FEL  oscillator  [2]  and  the  associated  wall-plug 
efficiency  with  an  ERL  [1,3].  The  optical  energy  is  PFEL  =  rjFEL  Ph  where  the  optical 

efficiency  in  the  low  gain  regime  is  tjfel  =  1/(2  Nw)  and  Nw  is  the  number  of  wiggler 
periods.  The  electron  beam  energy  exiting  the  decelerating  RF  fields  within  the  ERL  and 
entering  the  beam  dump  is  Pdump  =  Ih  (Emm  +  5Emc!  2) ,  where  SEmc  is  the  induced 
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electron  beam  energy  spread.  In  an  FEL  oscillator  the  energy  spread  is  proportional  to 
the  optical  energy,  a  Efel  ,  where  a  >  1 .  The  electron  beam  energy  spread  entering  the 
ERL  can  be  reduced  by  a  factor  ft  and  is  SEmc  =  (5 CX7Jfel  Eh .  The  phase  space  cooling 
parameter  /?  accounts  for  the  reduction  in  energy  spread  due  to  the  spatial  spreading  of 
the  beam  outside  of  the  wiggler.  The  wall-plug  efficiency  associated  with  the  oscillator 
is 


Vwp  ~  Gt]fel  — 


W  FEL 


(2) 


V FEL  (1  ■*"  &  2)  +  /L,  I  Ph  +  TJ  Pcryo  /  Pb 

As  an  illustration,  we  take  Nw  =  25  ,  a  =  4 ,  (3  =  1 ,  Eh  =  80  MeV ,  Ib  =  1  A , 

Pcryo  =  1  MW  ,  tj  =  0.6 ,  and  Emm  =  2  MeV  and  find  that  the  FEL  oscillator  interaction 
efficiency  is  tjfel  =  2%,  and  the  wall-plug  efficiency  is  ijWP  -  10%  .  Since  the 


oscillator  efficiency  is  ~  1/(2  Nw),  reducing  Nw  can  increase  the  wall-plug  efficiency. 
However,  in  the  low  gain  regime  the  gain/pass  is  proportional  to  Nl  and  therefore  Nw 
cannot  be  reduced  below  the  oscillator’s  start  current  condition. 

b)  Amplifier 

The  efficiency  of  an  FEL  amplifier  can  be  high  in  the  high-gain  regime  in  which 
the  radiation  initially  grows  exponentially.  The  efficiency  can  be  further  increased  by 
tapering  the  wiggler  or  by  frequency  detuning.  In  addition,  pulse  slippage  is  significantly 
reduced  compared  to  the  low  gain  oscillator  regime  [4],  Furthermore,  in  a  high-gain 
uniform  wiggler  amplifier  the  radiation  can  be  optically  guided  if  the  electron  current  is 
sufficiently  high  [4],  In  a  tapered  wiggler  amplifier,  optical  guiding  can  also  be  achieved 
even  in  the  trapped  particle  regime. 
i)  Uniform  Wiggler 

The  wall-plug  efficiency  associated  with  the  untapered  wiggler  amplifier  is 

n  = _ USm, -  (3) 

0  +  afim  +  C,  If, ;  +>lPa»IP> 

where  the  electron  beam  power  entering  dump  is  Pdump  =  lh  ( Emin  +  SEamp  /  2)  and  the 


induced  energy  spread  is  8Eamp  =  ParjFELEh.  The  FEL  interaction  efficiency,  rjFEL ,  in  a 
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uniform  wiggler  amplifier  can  be  increased  by  approximately  a  factor  of  ~  1.5  by 
frequency  detuning  [4], 

it)  Tapered  Wiggler 

The  wall-plug  efficiency  in  the  tapered  wiggler  amplifier  configuration  depends  in 
detail  on  the  bunch  dynamics  in  the  ERL.  This  analysis  is  presented  in  Section  III. 
However,  by  making  a  number  of  simplifying  assumptions  an  approximate  expression  for 
the  wall-plug  efficiency  can  be  obtained.  A  more  accurate  values  for  the  wall-plug  is 
presented  in  Section  IV. 

In  a  tapered  wiggler  amplifier  the  FEL  output  power  is  [4,  5] 

^FEL  ~  ^FEL^b  ~  ^Itrap  ^b  '  (4) 

where  Jj FEL  =  T]trap  APh/  Pb  is  the  FEL  interaction  efficiency,  rjlrap  is  the  trapping 

efficiency,  i.e.,  fraction  of  electrons  trapped  in  the  ponderomotive  bucket,  and 
A Ph  =  IhAEb  and  A Eh  is  the  change  in  the  electron  beam  energy  (tapered  wiggler) 

associated  with  the  trapped  electrons.  Figure  2  shows  a  sketch  of  the  electron  energy 
distribution  in  a  tapered  wiggler  amplifier  indicating  the  various  stages  of  the  trapped  an 
untrapped  electrons.  If  both  the  trapped  and  untrapped  electrons  are  equally  decelerated 
and  the  energy  spreads  are  unchanged  in  the  ERL,  the  electron  beam  energy  entering  the 
beam  dump  isEdump=Tjlrap(Emm+SE,mp/2)+(l-T]lrap)(AEh  +  Emin+SElrap/2),  which  can  be 

rewritten  as 

E dump  ~  Q  ~  V trap}  ^b  +  ^min  +  ^ trap  ^  ^  i  (5) 

where  SEtrap  is  the  induced  full  energy  spread  associated  with  the  trapped  electrons.  Note 

that  to  this  level  of  approximation  Edump  is  independent  of  SEunlrap .  The  energy  spread 
associated  with  the  trapped  electrons  is 

1/2 

a'J2(Eh-AE„),  (6) 

where  aR  =  6xlO“10  >l[pm]/^[W/cm2]  is  the  normalized  vector  potential  associated 
with  the  optical  field  of  wavelength  A ,  and  aK  is  the  wiggler  strength  parameter.  For  a 
high  average  power  FEL,  the  peak  optical  intensity  within  the  wiggler  is 
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1  peak  =  2 Ppeak  / (k R2l)  ~  1010  W/cm2  (aR  ~  10-4)  where  /?L  is  the  laser  spot  size,  and 
the  fractional  trapped  electron  energy  spread  is  SEtmp  I  Eh  ~  10”2 . 

The  wall-plug  efficiency  associated  with  the  tapered  wiggler  FEL  amplifier 
becomes 


Vwp 


_ Thsi _ 

H PEL  I  ?1  trap  ^min  +  ^trnp  ^  +  7  PcrytJ 


where  we  used  the  relation  in  Eq.(4). 
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III.  Dynamics  of  a  Charge  Bunch  in  an  ERL 

The  longitudinal  and  transverse  dynamics  of  a  relativistic  charge  bunch  are 
described  by  a  set  of  coupled  envelope  equations  for  the  root  mean  square  (rms)  bunch 
length  and  radius.  The  bunch  length  and  radius  dynamics  are  determined  by  the  effects 
of  the  RF  linac  field,  space  charge  fields,  inductive  (self  magnetic)  fields,  external 
focusing  fields  as  well  as  the  axial  and  transverse  energy  spreads  (emittances). 

The  space  charge  and  RF  fields  acting  on  the  electron  bunch  are  given  in  the 
Appendix.  The  axial  and  radial  self  electric  fields  with  respect  to  the  center  of  the  bunch 
are 

EScAz,»  =  Juf/t  M(DU  -  MO).  (8a) 

5V5  Kh  Lh 

E^r’,)  =  iht^rym)]r'  (8b) 

where  the  function  M(£)  depends  on  the  bunch’s  aspect  ratio,  energy  and  surrounding 
geometry  (see  Fig.  3  and  the  Appendix),  q  N  is  the  bunch  charge,  Lh  and  Rb  are  the  rms 


bunch  length  and  radius  in  the  laboratory  frame  and  zc(t)  is  the  axial  position  of  the 


bunch  centroid.  In  the  absence  of  walls  M(^)  is  a  function  of  the  aspect  ratio  and  energy 


only. 


M(<f)  = 


1-f 


( 


2<f 


V 


In 

v 


(9) 
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where  %  =  (1  —  Rb  l(yLbj1)'12 .  The  express  for  M(^)  in  Eq.  (9)  is  valid  for  the  full 
range  of  Rh  l(yLb) .  However,  for  the  parameters  considered  in  this  paper 
Rh/(yLb)«\,  Rh~ Lh,  and  M(<f)  -  -f  (0.31  +  ln£). 

The  ERL  fields  are  given  by 

ERFtZ(z, 0  =  Erf  sin(^f),  (10a) 

ERFr(r,z,t)  =  -ERF(kRFr/  2)cos(4/f),  (10b) 

BRFQ(r,z,t)  =  -ERF((0RFr/ 2c)cosC¥c) ,  (10c) 

where  £R/,  is  the  field  amplitude,  4*c  =  (z')  -  0)RF  lePc(z'))dz'  +  4^(0)  is  the 

o 

phase,  coRF  is  the  RF  frequency  and  kRF  =  2 Kl  XRF  is  the  local  wavenumber  of  the  RF 
field. 

The  envelope  equations  are  derived  in  the  Appendix  and  are  given  by 


d% 

dz 2 

d% 

dz2 


( 


1 


\ 


v3+  PlfcjYc  dz  dz 


1  d/c  r)Lh 


Kl  A  ~ 


Pi  t  E\ 


=  0, 


-  K\Rb 


A  Yc  dz  3 z 


Alt  Rt 


=  0, 


(11a) 

(lib) 


where  K2  =  K2RF  cos  4^  +  K^z,  K\  =  K2cS  -  K2  -  K2RF^ 


K2cz=3/(5y[5)q2N/(mc2)M^)/(j32fcR2hLb),  K2RF  =  qERFkRF /(me2  p2  fc), 

Kl.x  =  3/(10V5)^2  N/(mc2)(l  -  M(<f;))/(p2  y]  R2b  Lh) .  The  quadrupole  focusing  term  is 
Kq  =  ^q  Bn/  rQ^/(yi Pc  me) ,  where  Bn  is  the  magnetic  field  and  rQ  is  the  characteristic 
gradient  scale  length  and  Bo  /  rQ  ~  5kG/cm  .  The  RF  focusing/defocusing  term  is 

kIf, i  =xqERF  /(me2  fc  p]  ARF)cos  4/f .  The  energy  associated  with  the  bunch  centroid 
is  given  by 


dyc  _  qEi 
dz 


RF 

me 2 


Ac*  in(4/J, 


(12) 


where  y  =  (1  -  p2) 


2  v  — l  /  2 
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IV.  FEL  Amplifier  Wall-Plug  Efficiency 

In  this  section  we  calculate  the  FEL  wall-plug  efficiency  for  a  tapered  amplifier 
configuration.  We  numerically  solve  Eqs.  (1 1)  and  (12)  using  electron  beam  parameters 
appropriate  to  a  high-power  FEL.  Figure  4  shows  a  tapered  FEL  amplifier  configuration 
in  which  the  electron  beam  from  an  RF  gun  is  accelerated  in  the  linac,  enters  a  tapered 
FEL  wiggler,  generates  FEL  radiation,  and  emerges  as  a  spent  beam  consisting  of  trapped 
and  un-trapped  particles.  The  spent  electron  bunch  is  shown  schematically  in  Fig.  5. 
Figure  5(a)  shows  the  initial  electron  beam  distribution  exiting  the  FEL  wiggler. 

Injecting  these  two  populations  at  the  same  phase  of  the  RF  wave  causes  them  to  be  de- 
accelerated  to  the  final  configuration  shown  in  Fig.  5(b)  in  which  the  trapped  particles 
have  a  minimum  energy  of  £%,  and  the  un-trapped  particles  have  a  minimum  energy  of 
~  £min  +  'HpEL^b  ■  The  un-trapped  particles  can  be  de-accelerated  to  £min  by  injecting 

them  at  a  different  phase  relative  to  the  trapped  particles,  i.e.,  one  for  which  the  de- 
accelerating  field  is  larger,  as  shown  in  Fig.  5(c).  This  increases  the  recovered  energy  by 
an  additional  amount  ~  TjFELEh,  which  increases  the  overall  wall-plug  efficiency. 

In  the  example  that  follows,  the  RF  linac  field  is  assumed  to  have  a  wavelength  of 
40  cm  and  an  accelerating  gradient  of  1 5  MeV/m.  The  accelerating/de-accelerating 
cavity  is  5.5  m  long.  An  external  quadrupole  focusing  field  with  normalized  amplitude 
qBo  /(a me2)  =  0.0 2 cm-2 is  assumed. 

The  initial  electron  bunch  before  acceleration  in  the  linac  is  assumed  to  have  the 
parameters  listed  in  Table  1,  e.g.,  energy  of  5  MeV  and  a  bunch  charge  of  1  nC.  Injecting 
the  bunch  at  an  initial  phase  of  =  1 .9 ,  it  is  accelerated  to  ~  83  MeV  at  the  end  of  linac 

as  shown  in  Fig.  6(a).  The  energy  spread  decreases  from  0.12%  to  0.008%  during 
acceleration  (Fig.  6(b))  and  the  bunch  length  remains  relatively  contstant  (Fig.  7(a)).  The 
bunch  radius  undergoes  oscillation  due  to  the  RF  and  external  focusing  fields  (Fig  7(b). 

Before  entering  the  wiggler,  the  beam  is  focused  to  a  spot  size  of  0.2  mm  and 
compressed  to  ~1  psec  to  increase  the  peak  current  for  the  FEL  interaction  in  a  tapered 
wiggler.  At  the  wiggler  exit,  the  bunch  length  is  ~  1  psec  and  the  trapped  and  un-trapped 
particles  are  assumed  to  have  a  longitudinal  emittance  ofe^  z  ~  2.5  MeV  -  psec  ~  2500 
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keV-psec.  This  corresponds  to  a  fractional  energy  spread  of  ~  3%.  In  a  properly  tapered 
wiggler,  ~  80%  of  the  particles  can  be  trapped. 

Before  re-entering  the  linac  to  be  de-accelerated,  it  is  assumed  that  the  bunch  is 
lengthened  to  5  psec  to  reduce  the  energy  spread.  The  trapped  particles  are  injected  at  a 
phase  4^  =  4. 1 8  and  the  un-trapped  particles  are  injected  at  4/0  =  4.73 .  This  causes  both 

populations  to  be  de-accelerated  to  a  final  energy  <  2  MeV  at  the  end  of  the  linac,  as 
shown  in  Fig.  6(a).  Figure  6(b)  shows  that  the  fractional  energy  spread  of  the  trapped 
particles  increases  to  >  30%  near  the  end  of  the  linac  before  decreasing  to  a  value  of  ~ 

7%  at  the  exit.  The  un-trapped  particles  have  a  final  energy  spread  of  ~  1 8%.  The  bunch 
length  of  the  trapped  and  un-trapped  particles  remains  relatively  constant  until  the  end  of 
the  linac  where  the  trapped  particles  cause  the  bunch  length  to  increase  to  ~  2.5%  of  the 
RF  wavelength  [Fig.  7(a)].  The  final  bunch  radius  of  the  trapped  particles  is  ~0.1  cm. 

Figure  8  shows  the  wall-plug  efficiency  for  oscillator  (green  curve),  un-tapered 
amplifier  (blue  curve),  and  tapered  amplifier  (red  curves)  configurations  versus  minimum 
beam  energy.  For  the  tapered  amplifier,  the  solid  red  curve  denotes  the  case  where 
trapped  and  un-trapped  particles  are  injected  at  the  same  phase,  *F0  =  4. 1 8 .  The  dashed 
curve  denotes  the  case  for  which  trapped  and  un-trapped  particles  are  injected  at  different 
phases,  i.e.,  4/0  =  4.18  for  trapped  particles  and  =  4.73  for  un-trapped  particles.  The 

parameters  for  Fig.  8  are  Eh  =  80  MeV ,  Ih  =  1  A ,  JJRF  =  0.6 ,  Pcryo  =  1  MW  ,  77  =  0.6 , 

£  =  0.3 ,  77inj  =  0.6 ,  and  /?  =  1 . 


V.  Conclusions 

We  have  calculated  the  wall-plug  efficiency  for  an  FEL  oscillator,  amplifier,  and 
tapered  amplifier  configuration  employing  an  energy  recovery  linac.  The  wall-plug 
efficiency  depends  on  various  FEL-ERL  parameters  and  can  be  improved,  for  example, 
by  reducing  the  minimum  energy  entering  the  beam  dump,  Emin .  By  reducing  Emin  from 
4 MeV  to  2  MeV  and  keeping  all  other  parameters  the  same  the  wall-plug  efficiency  for 
the  oscillator  increases  from  10%  to  12%;  for  the  uniform  wiggler  amplifier  the 
efficiency  increases  from  10%  to  14%;  and  for  the  tapered  wiggler  amplifier  the 
efficiency  increases  from  22%  to  24%.  For  typical  MW  class  FEL  parameters,  we  find 
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that  the  wall-plug  efficiency  for  the  oscillator  is  in  the  range  tjwp  ~  1 0  - 1 5  %  For  the 
uniform  wiggler  amplifier  tjwp  ~  10-15%,  and  for  the  tapered  wiggler  amplifier 
J?WP  *  20-25% 

Acknowledgments:  The  work  was  sponsored  by  the  Office  of  Naval  Research  and  the 
Joint  Technology  Office. 

Appendix 

The  longitudinal  and  transverse  dynamics  of  a  relativistic  charge  bunch  are 
described  by  a  set  of  coupled  envelope  equations  for  the  root  mean  square  (rms)  bunch 
length  and  radius.  The  bunch  length  and  radius  dynamics  are  determined  by  the  effects 
of  the  RF  linac  field,  space  charge  fields,  inductive  (self  magnetic)  fields,  external 
focusing  fields  as  well  as  the  axial  and  transverse  energy  spreads  (emittances). 

I.  Fields 

a)  Space  Charge  Fields  [6-8] 

The  potential  associated  with  a  stationary  (beam  frame)  elliptical  charge  bunch  of 
uniform  density,  with  radial  and  axial  semi-axis  of  (a,  zm  )  in  the  beam  frame,  centered  at 

Z  =  Ois 

0(r,z)  =  -Inpi—  —~^Hr2  +  M(£)zO,  (Al) 

2  J 

where  the  charge  density  is  po  =  (3/4 7r)qN /(a2  zm) ,  q N  is  the  charge  in  the  bunch. 

The  uniform  charge  ellipse  is  defined  by  r 2  la2  +  z2  /  z2m<  1 .  The  electro- static  space 
charge  field  in  the  beam  frame  is  EJC  =  -  VO .  In  the  beam  frame  the  function  M  (£) 
depends  on  the  aspect  ratio  and  geometry.  In  the  absence  of  walls  M  (£)  is  a  function  of 
only  the  aspect  ratio,  g  =  (1  -  o2  / z2m)U2 . 

In  the  laboratory  frame  the  space  charge  fields  are  modified  by  both  relativistic 
and  inductive  effects.  In  the  laboratory  frame  where  the  bunch  is  moving  with  constant 
axial  velocity  v ,  the  scalar  potential  is,  in  the  Lorentz  gauge,  given  by 
(v2  +  (\/ y1)d2 /d?j2)^(r,T])  =  - 4zrp(r,jj )  where  7]  =  z-vt ,  ^=(1-  v2/c2)“l/2  and 
p  is  the  charge  density.  In  the  laboratory  frame  there  exists  an  axial  vector  potential 
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given  by  A  =  (v/c)<J>  which  results  in  an  inductive  axial  electric  field.  The  rms  bunch 


radius  and  length,  in  the  laboratory  frame,  are  given  by  Rh  =  l^r2^  =  a/^5  and 

Lh  =  =  zm  /(V 5  y) ,  respectively.  The  scalar  potential  in  the  laboratory  frame  is 

found  from  the  expression  for  the  potential  in  the  beam  frame,  Eq.  (Al),  by  replacing  the 
variable  z  with  yrj  and  the  density  with  p0  =  (3/47r)(l/lO^)qN /(Rb  Lb) .  The 
potential  in  the  laboratory  frame  becomes, 

™ 

where 

1_A 2f  i  +  A  ] 

M(£)  =  —fr  In  — f  -  24  ,  (A3) 

v  v1  _  Cy  J 

and  %  =  (1  -  ( Rh  /  yLh)2)112 .  The  expression  for  M(£)  is  valid  for  the  full  range  of  the 
ratio  Rh  /(y  Lh )  and  is  shown  in  Fig.  3.  The  axial  and  radial  self  electric  fields  of  the 
bunch,  centered  at  r  -  0,  z  =  zt  -  vc  t ,  in  the  laboratory  frame 


(A4a) 


(A4b) 


b)  ERL  Fields 


The  ERL  electric  and  magnetic  fields  near  the  axis  are  given  by 


ERFJz,t)  =  Erf  sinl  \kRF(z)dz  +  <nRFt  , 

ERF,r(z,t )  =  -ERF(kRF  r / 2)cos(  jkRF(z)dz  +  (DRFt  , 


(A5a) 


(A5b) 


F  Z 

BRFS(r,z,t)  =  -ERF(G)RFr/2c)cos  jkRF(z)dz' +  (£>RFt  ,  (A5c) 


Erf  is  the  field  amplitude,  0)RF  is  the  RF  frequency  and  kRF  =  2n I XRF  is  the  local 
wavenumber  of  the  RF  field. 


II.  Bunch  Dynamics 
a)  Longitudinal  Bunch  Dynamics 

The  relativistic  longitudinal  electron  motion  is  given  by 


a2 


dt2 


f 


m 


fz  \ 

\ 

Erf  sin 

\dz  ( z )  “  (Orf  t 

+  Escjz,t) 

V  ) 

J 

(A6) 


where  z{t )  is  the  axial  electron  position  and  y(t)  =  (l  -  ( dz/dt )2  / c 2)  1/2 .  The  axial 
electron  position  can  be  written  as  z(t)  =  zc(t )  +  S z(t) ,  where  Sz  is  the  axial  position 
relative  to  the  centroid.  The  equation  of  motion  of  the  centroid  of  the  bunch  is  given  by 


a2zt.(0  _  q 

dt2 


RF 


sin^CO), 


(A7) 


c RF 


while  Sz{t)  is  given  by 

d2Sz_o;if_ 


dt 2 


VRF 


/ 

r  > 

\ 

sin 

+  \dzkRF(z’) 

-sin(*Fc) 

V 

l  *  J 

J 

±llsc,z°z  T"  TT’ 

Yc  ot  at 


(A8) 


where  '¥c(t)=  \dzkRF{z)  -  0)RFt  +  ¥c(0) ,  Q.2SCZ  =  3/(5  J5)(q2  N  /  m)M  !{fc  R2hLh ) , 
0 


=  <7  £rf  ^rf  m)  >  ^RF  =  ^RF  <X(0) ,  7,  =  0  -  A2)  1/2  and  A  =  c  1 / at .  The 
energy  associated  with  the  bunch  centroid  is  given  by 

^  jBc  sin^.)  =  ^—S£-Pcs\ni^Vc).  (A9) 

at  ckRf  mc 

For  kRF  constant  and  \kRF  Sz\  «  1  we  can  write 


d2Sz 

dt2 


kl]5z 


,  3  dydSz  =  () 
yc  dt  dt 


(A10) 


where  Q2  =  Q2RF  cos T,.  +  Q2c, ,  Tc  =  kRFzc  -0)RFt,  zc  =  qERF  Fyfc  m)sin(vf/t  )  and  the 


overdot  indicates  d/dt . 

The  longitudinal  electron  motion  about  the  centroid  is  given  by  Eq.  (A  10),  which 
upon  multiplying  by  Sz  and  Sz  and  combining  the  resulting  averaged  equations  leads  to 
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3^  +  ^)(44  +  (4)2  -&A  +  ir: thu) - &UK  =  o,  (aid 


where  L,,  =  (^5  z2}  is  the  rms  bunch  length.  It  can  be  shown  that  the  above  equation 

can  be  written  as  a  total  derivative, 

3 


-Q?3  +  3/;7,44))  =  o. 


(A  12) 


and  t  U,  ~  &  Ll  +  3  f'  ye  4  4)  =  c 2  ,  where  en_z  is  the  constant  normalized 

longitudinal  emittance.  Since  LhLh  -  D.I  +  3  y~l  yc  =  (^Sv2^  —  ,  we  can  write 

the  emittance  term  as  c2  £2 ,  =  yf  L^i^Sv2^  -  X)2)  where  Sv  =  Si.  Assuming  self¬ 
similar  behavior,  Sv  =  SzLh/Lh  +  vth ,  the  normalized  longitudinal  emittance  becomes 

el=  =  i^((*'!)-(Z‘)Z)=^^4  =  £-r-?i-  <AI3> 


C  C 

The  longitudinal  envelope  equation  for  the  bunch  is 

^  4  +  —  _ 


c2  f„2 . 


=  o. 


(A  14) 


dt2  ~~z  yc  dt  dt  tfL\ 

The  normalized  longitudinal  emittance  en  .  =  Lh  ythl  fic  is  a  Lorentz  invariant  which  is 

also  commonly  written  as  enz  =e*nz/mc  where  (Vmc)  ~  5.87x1 04[cm/(eV-sec)], 

£* .  =  A Eh  Arh  [eV-sec]  and  Aih  =  Lb/vc  is  the  bunch  duration  [9]. 

b)  Transverse  Bunch  Dynamics 

Using  a  similar  approach  an  envelop  equation  for  the  rms  bunch  radius  can  be 
obtained.  The  radial  electron  motion  about  the  center  of  the  bunch  is  given  by 
Sr  -  Q2±Sr  +  y~x  yc8r-  0 ,  where  Q2  =  Q2  x  -  Q2 , 

£22f±  =  3/(10a/5)^2  N /(y^  m)(\  -  Lh)  is  due  to  the  transverse  self  fields  and  Q2 

represents  external  transverse  defocusing  or  focusing  due  to  the  RF  and  quadrupole 
fields.  The  transverse  envelope  equation  is 


aX  1  dy  d_R, 
dt2  yc  dt  dt 


h  -  niR„ 


2  2 

i3 


(A  15) 
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where  enL  is  the  rms  normalized  transverse  emittance. 

c)  Change  of  Variables 

It  is  useful  to  employ  the  axial  coordinate  z  as  the  independent  variable  instead 
of  the  time  t .  Performing  this  transformation  the  equations  become 


,  2  i  1  1  1  *Yc 

dz2  l  PcfcjYcdzdz 


K]  L. 

z  b 


(A  16a) 


,  _i_an 

dz2  p2ycdzdz 


d  Yc  d  Eftf  n 

~\  1  r'c 

dz  me 


sin  ('Pc), 


(A  16b) 
(A  16c) 


where  K2  =  K2RF  cos^  +  K^z,  K*  =  K2^  -  K 2R?1  -  K2Q, 
K2CtZ  =  3/(5  V5 )(q2  N/m)M/(c 2  ft]  fc  R2Lb ) , 


Z 

%  =  ]{kRF(z)  ~  coRF  l(cpc(z)))dz  +  T(0) , 

0 

K2Ct±  =  3/(10V5)^2  N/(mc2){\  -  M  (£))/(#  fc  R2  4) ,  K2RF  =  qERF  kRF  /(me2  p2  fc). 


kIf,l  =  zrqERF  /(me2  fc  $  ARF )eos'Fc  is  due  to  RF  focusing/defocusing,  and 

Kq  =  \qBa/rQ\/(yc Pcme)  is  due  to  quadrupole  focusing,  Ba  is  the  magnetic  field  and  rQ 

is  the  characteristic  gradient  scale  length. 
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Bunch  parameters 

Initial 

bunch 

Trapped  particles 
after  wiggler 

Un-trapped  particles 
after  wiggler 

Energy  (MeV) 

5 

70 

83 

Longitudinal  emittance 
[keV-psec] 

100 

2500 

2500 

Transverse  emittance 
[mm-mrad] 

5 

5 

5 

Bunch  length  [psec] 

15 

5 

5 

Energy  spread  [%] 

0.12 

0.7 

0.6 

Radius  [mm] 

1 

0.2 

0.2 

Bunch  charge  [nC] 

1 

0.8 

0.2 

Injection  phase  [rad] 

1.9 

4.18 

4.73 

Table  1 :  Electron  bunch  input  parameters  used  in  the  numerical  examples  of  Figures  6 
and  7. 
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Figure  1 :  Power  flow  diagram  of  an  FEL  employing  an  ERL  which  is  used 
to  derive  the  wall-plug  efficiency. 
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Figure  2:  Illustration  of  the  electron  energy  distribution  in  a  tapered  FEL  amplifier,  (a) 
shows  the  electron  distribution  function  at  the  wiggler  entrance  ( z  =  0 ),  a  beam  of  N 
electrons  with  energy  Eh .  (b)  shows  the  electron  distribution  at  the  tapered  wiggler  exit 
(z  =  Lw).  A  fraction  of  electrons  rjlrap  are  trapped  and  de-accelerated  to  an  average 
energy  Eh  -  A Eh .  The  trapped  particles  have  an  induced  energy  spread  <%,rap .  (c)  shows 

the  electron  distribution  at  the  beam  dump.  Both  trapped  and  un-trapped  electrons  are  de- 
accelerated  by  an  amount  -  A Eh  -  (SElrap/2)  -  E^  .  The  total  energy  entering  the 

beam  dump  in  (c)  is  given  by  Eq.  (5) 
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Figure  3:  Plot  of  the  space  charge  function  M  versus  Rh/(yLh),  Eq.  (9).  The  self-fields 
of  an  ellipsoidal  charge  with  length  Lh  and  radius  Rh  are  expressible  in  terms  of  M  as 
given  in  Eq.  (8).  In  (a)  the  function  M  is  plotted  for  0  <Rh RyLh)  <  1 ,  while  in  (b),  the 
range  of  Rh  RyLh)  is  extended  to  20. 
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Figure  4:  Schematic  diagram  showing  the  electron  beam  from  an  RF  gun  being 
accelerated  in  the  linac,  entering  a  tapered  FEL  wiggler,  generating  FEL  radiation,  and 
emerging  as  a  spent  beam  consisting  of  trapped  and  un-trapped  particles.  The  trapped 
and  un-trapped  particles  are  then  injected  into  the  linac  at  different  phases  and  de- 
accelerated  to  recover  most  of  their  energy  before  being  sent  to  a  beam  dump. 
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Figure  5:  Schematic  diagram  showing  the  trapped  (red)  and  untrapped  (blue)  electron 
energy  distributions  (a)  at  the  end  of  the  tapered  wiggler,  and  after  energy  recovery  when 
(b)  the  particles  are  injected  at  the  same  RF  phase,  and  (c)  the  trapped  and  untrapped 
particles  are  injected  at  different  phases  in  order  to  de-accelerate  them  to  the  same 
minimum  energy  and  enhance  wall-plug  efficiency. 
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Figure  6:  Beam  energy  (a)  and  energy  spread  (b)  versus  longitudinal  position,  z,  within 
the  linac.  Curves  denote  the  electron  bunch  before  the  FEL  interaction  (black  curves), 
and  the  trapped  (solid  red  curves)  and  untrapped  (dashed  red  curves)  particles  after  the 
FEL  interaction  in  a  tapered  amplifier.  Parameters  are  listed  in  Table  1 . 
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Figure  7:  Bunch  length  normalized  to  RF  wavelength  (a)  and  bunch  radius  (b)  versus 
longitudinal  position,  z,  within  the  linac.  Curves  denote  the  electron  bunch  before  the 
FEL  interaction  (black  curves),  and  the  trapped  (solid  red  curves)  and  untrapped  (dashed 
red  curves)  particles  after  the  FEL  interaction  in  a  tapered  amplifier.  Parameters  are  listed 
in  Table  1. 


23 


Figure  8:  Wall-plug  efficiency  for  oscillator  (green  curve),  un-tapered  amplifier  (blue 
curve),  and  tapered  amplifier  (red  curves)  configurations  versus  minimum  beam  energy. 
For  the  tapered  amplifier,  the  solid  red  curve  denotes  the  case  where  trapped  and  un¬ 
trapped  particles  are  injected  at  the  same  phase,  'Fq  =  4. 1 8 .  The  dashed  curve  denotes  the 
case  for  which  trapped  and  un-trapped  particles  are  injected  at  different  phases,  i.e., 

=  4. 1 8  for  trapped  particles  and  M',,  =  4.73  for  un-trapped  particles.  ERL  parameters 

are  Eh  =  80  MeV ,  Ih  =  1  A  ,  tjrf  =  0.6 ,  Pctyo  =  1  MW  ,  rj  =  0.6 ,  £  =  0.3 ,  T]mi  =  0.6 ,  and 
0  =  \. 
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